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Using a scalar field as an intrinsic ’clock’ while investigating the dynamics of gravitational systems 
has been successfully pursued in various researches on the border between classical and quantum 
gravity. The objective of our research was to check explicitly whether the scalar field can serve 
as a time variable during dynamical evolution of the matter-geometry system, especially in regions 
of high curvature, which are essential from the perspective of quantum gravity. For this purpose, 
we analyzed a gravitational collapse of a self-interacting scalar field within the framework of general 
relativity. The obtained results indicated that the hypersurfaces of constant scalar field are spacelike 
in dynamical regions nearby the singularities formed during the investigated process. The scalar 
field values change monotonically in the areas, in which the constancy hypersurfaces are spacelike. 


I. INTRODUCTION 

Canonical gravity is an approach to general relativity 
based on an arbitrary slicing of spacetime into spacelike 
hypersurfaces labeled by a single parameter, which can 
be interpreted as a time variable. All possible slicings are 
physically equivalent. The dynamical evolution in such 
a system is governed by the Hamiltonian, which is also 
a constraint and a generator of gauge transformations 
between the spacelike slices. The quantizations of grav¬ 
ity based on its canonical formulation encounter the so- 
called problem of time. It is related to the fact that there 
is no clear notion of time in gravitational systems, which 
could be transferred from the classical to quantum level. 
The problem reveals itself as crucial in attempts to define 
the dynamics of a quantized system. 

The idea of using matter to quantify time during 
the course of gravitational evolution was proposed by De- 
Witt in [1] . In the context of quantum theory of gravity, 
it was suggested that the matter content of a spacetime 
can play a role of a ’clock’ when investigating the dynam¬ 
ical behavior of geometry. One of the degrees of free¬ 
dom of the system is interpreted as a dynamical ob¬ 
server. It serves as a reference for the remaining degrees 
of freedom, whose dynamics is followed with respect to it. 
Potentially, this concept can be also applied inversely, 
i.e., the geometry could be used to introduce a physical 
time when the dynamical evolution of matter is exam¬ 
ined. 

The notion of internal time is especially important 
in non-perturbative quantum gravity, which deals with 
a problem of investigating time evolution without fixed 
geometry of a background spacetime [2]. A Hamilto- 
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nian which describes quantum evolution of a gravita¬ 
tional field with respect to time measured by a ’physi¬ 
cal clock’ provided by spacelike constancy hypersurfaces 
of a scalar field was presented in m- It was defined 
through a regularization procedure based on the Ashtekar 
variables HE] and the loop representation laiT] and was 
shown to be finite and diffeomorphism invariant. Re¬ 
cently, the idea was developed by applying loop quan¬ 
tum gravity tools. The scalar field, which is an intrin¬ 
sic element of the model, was treated as the time vari¬ 
able for the relational Dirac observables, whose dynamics 
could be traced with respect to the field ElEI. 

The scalar field was also used to quantify the pas¬ 
sage of time during quantum cosmological evolutions. 
The studies of the inflation epoch of the Universe 
in the regime of loop quantum gravity was studied with 
the use of a scalar field with an arbitrary potential HQ). 
The gauge for the Hamiltonian constraints was chosen 
so that the scalar field was constant on spacelike hyper¬ 
surfaces. Treating the scalar field as a time variable al¬ 
lowed writing the Hamiltonian constraints in the form 
of a time-dependent Schrodinger equation. A similar ap¬ 
proach and a massless scalar field is used in the loop 
quantum cosmology munj. 

The Brans-Dicke scalar field, which is non-minimally 
coupled to geometry, was also examined as a ’clock’ 
in the quantum cosmological context m- It was claimed 
that since the field is a monotonic function of cosmo¬ 
logical time, it may serve as an internal time variable. 
A spatially fiat and isotropic cosmological Brans-Dicke 
model was quantized within the framework of loop quan¬ 
tum cosmology and the effective Hamiltonian was thus 
obtained. 

A homogeneous, massless, minimally coupled scalar 
field was used to calculate the tunneling decay rate m 
of a simple harmonic universe m- It is a cosmological 
model, which oscillates due to the presence of a negative 
cosmological constant and a matter component with re- 
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pulsive gravity. An additional field with negligible contri¬ 
bution to the total energy density of the system, which 
served as a ’clock’, was included within it min]. Its 
presence was necessary, because the scale factor, which 
is the only dynamical variable of the model, is not mono¬ 
tonic in an oscillating universe and hence it cannot be 
treated as a time variable. 


The version of the Wheeler-DeWitt equation, which 
uses the scalar field as a variable common for internal 
time and a Hamiltonian evolution parameter, and its 
quantum counterpart for the timelike case were presented 
in m- The former describes the evolution of four¬ 
dimensional Lorentz hypersurfaces embedded in a five¬ 
dimensional spacetime, along the massless scalar field. 
The latter has the form of the Schrodinger equation, 
i.e., it contains the internal time derivative on the right- 
hand side and presents the evolution with respect to it. 

Apart from the scalar field, dust and radiative fluid 
were also used to provide a matter degree of freedom, 
which allowed tracing the temporal evolution of a gravi¬ 
tational system. The comoving coordinates of dust parti¬ 
cles and the proper time along the dust world lines served 
as canonical coordinates in the phase space of the dy¬ 
namical system [19]. The issue of unitarity of the evolu¬ 
tion proceeding in the isotropic and homogeneous Brans- 
Dicke quantum cosmological model, with a time variable 
defined by the matter fluid, was discussed in [20] . 

Despite the important achievements in formulating 
a consistent description of a time evolution of a gravi¬ 
tational system, which were enabled by treating matter 
as a time variable, a question about the relevance of such 
an approach remains open. In general, it is assumed 
that the behavior of matter is desirable during the inves¬ 
tigated part of an evolution, i.e., that the spacelike char¬ 
acter of the chosen spacetime slices is preserved and that 
the selected time parametrization remains monotonic 
during its course. However, the arguments in favor 
of these assumptions are limited to certain cases (such 
as homogeneity of a scalar field during an initial phase 
of inflation, which results in a spacelike character of con¬ 
stant field hyper surfaces m) and there does not exist 
a general justification for them. While for perturbations 
of the cosmological spatially homogeneous solutions it is 
plausible that the constancy surfaces of the scalar field 
are still spacelike, a character of those surfaces is highly 
non-trivial in the case of the spherically symmetric gravi¬ 
tational collapse. The latter is also a hot subject of quan¬ 
tum gravity [21]. 


In the present paper we considered an evolution 
of a scalar field, which was its dynamical gravitational 
collapse, within the framework of a four-dimensional 
formulation of general relativity. The main objective 
of the research was to examine the feasibility of using 
the scalar field as a time variable in spacetime regions 
of high curvature, which are essential for the applications 
of the quantum theory of gravity. The problem was posed 
and its theoretical formulation was introduced in Sec¬ 
tion |TT] Section [ml contains the description of numerical 


algorithm used during computations, the specifics of its 
implementation and its tests. The results were presented 
and discussed in Section [TV] while conclusions were gath¬ 
ered in Section IVl 


II. THEORETICAL FRAMEWORK 

The Lagrangian of the scalar field ^ is 

£ = ( 1 ) 

and its stress-energy tensor has the form 

T^u = ( 2 ) 

where is a metric tensor of spacetime. The equation 
of motion for the scalar field derived from the variational 
principle is 


= 0. (3) 

The gravitational collapse was investigated in spherical 
symmetry. Since the evolving field is massless, it is con¬ 
venient to follow the process in double null coordinates 
= {u^v^O^cj))^ in which the line element has the form 

ds^ = —a^dudv + r‘^dQ?, (4) 

where a and r are arbitrary functions of the retarded time 
u and the advanced time while dCl‘^ = d0‘^ + sin^ Odcj)^ 
is the line element of the unit sphere, 0 and (p are an¬ 
gular coordinates. The employed method is referred 
to as the 2-T2 formalism [22] . 

The scalar field equation ^ written in the selected 
coordinates is 


r^,uv + = 0 (5) 


and the Einstein equations for the gravitational field are 
the following: 


a 


MV T ^ M^M T . — O5 


r 


M 

.aM,u , ^^2 




(6) 

( 7 ) 

(8) 
( 9 ) 


where denotes the partial derivative. 

The system of equations (§ § describes 

the spherically symmetric dynamical gravitational col¬ 
lapse of a scalar field in double null coordinates. The par¬ 
ticulars of the numerical algorithm employed to solve 
the above equations and its tests are presented in Sec¬ 
tion [TTIl 

The only arbitrary initial data of the computations is 
the profile of the evolving scalar field, which is imposed 
on the null u = const, hypersurface. The initial hyper¬ 
surface will be denoted as = 0 for further convenience. 
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The investigated families of initial profiles are presented 
in Fig. They are one-parameter classes of functions 
parametrized by the amplitude A^. The values of 
used in the plots were chosen so that the collapse resulted 
in a singular spacetime, as will be explained in more de¬ 
tail below. From now on, the profiles presented in Fig.[^ 
with the amplitudes specified in its caption will be re¬ 
ferred to as profiles PI to P9. The profiles describe spher¬ 
ically symmetric shells of matter as they tend to zero 
for ^ oc and —> 0. The former condition ensures 

asymptotic flatness of spacetime. The latter is equiv¬ 
alent to vanishing as r ^ 0, because the gauge free¬ 
dom of equations in i;-coordinate was thus fixed, see Sec¬ 
tion HIl 


III. NUMERICAL COMPUTATIONS 


The considered dynamical gravitational collapse is de¬ 
scribed by second-order differential equations 
In order to solve the system numerically a set of auxiliary 
variables was introduced 


c = 


a 


J ^5^ 

d = 

a 


f = 

X = 


9 = 

y = 


( 10 ) 


The quantities / and g are related to the expansions Ou 
and Oy of the null geodesic vector fields du and, respec¬ 
tively, dy via Ou = 2r“^/ and Oy = 2r~^g [23]. The hy¬ 
persurfaces determined by / = 0 and ^ = 0 are the trap¬ 
ping and anti-trapping horizons in the spacetime, respec¬ 
tively [24]. The squares of the quantities x and y repre¬ 
sent fluxes of the scalar field through null hypersurfaces 
of constant v and respectively. The substitutions (10) 
allowed us to obtain the system of first-order differential 
equations, which determines the course of the investi¬ 
gated process 


U1 : — ac = 0, (11) 

V2 : a^y - ad = 0, (12) 

V3:r,u-f = 0, (13) 

U4:r,-^ = 0, (14) 

U5:4>,^-x = 0, (15) 

-^ = 0, (16) 

El : f,u = ‘^cf - (17) 

-B2 ; = "Idg - ry^, (18) 

E3 : = -Ar“\ (19) 

E4 : = d^u = Ar“^ — xy, (20) 

S : rx^y = ry^u = -V, (21) 
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where A = ^ -h and rj = gx fy. The evolution of d 
and y along u was governed by relations EA and S', re¬ 
spectively. The remaining quantities, that is a, 4>, r, / 
and X, evolved along the ^’-coordinate according to equa¬ 
tions U2, U6, E2^ U4, E3 and S. The function c did not 


play an active role during the computations, so it was 
ignored. 

The evolution was followed on a two-dimensional grid 
constructed in the (r’i4)-plane with integration steps in v 
and u directions denoted as hy and hu^ respectively, 
and initially equal to 10“^. The numerical algorithm 
proposed in [25] was implemented. The general form 
of equations (11) - © can be written as 

f « = F (f, g), g,„ = G (f, g). (22) 

Calculating a value of the particular function f {u, v) 
or g{u^v) at the point {vq^uq) required the knowledge 
of the values of functions F and G at points (i;o — hy^uo) 
and {vq.uo - hy) 


f| 




(vo,uo) 


+ f| 


{VQ,UQ—hu) 


+ 


+/i„F(ff,gg) 




(vo,wo) ^ ^ 

-\-hyG (ff. 


{vQ-hy.uo) 


{vO,Uq) 


(23) 


(24) 


where two types of auxiliary quantities were introduced 


ffl 


gg 


(vo,wo) 

+h^F(f, g) 




+ 


(^^0,Wo) 




I {vo,UQ-hu)^ 

{vo — hv,uo) 


-y(G(f,g) 




G(ff,gg) 


(^^0,Wo) 


(25) 


(26) 


Double null coordinates ensured regular behavior of all 
the evolving quantities within the domain of integration 
except the vicinity of r = 0. During numerical computa¬ 
tions considerable difficulties also arose nearby the event 
horizon, where the function / was discontinuous. A dense 
numerical grid was necessary to conduct the calculations 
in the neighborhood of the event horizon and to examine 
the behavior of fields beyond it, especially for large val¬ 
ues of advanced time. Since regions in which the dense 
grid was required were well-defined within the domain 
of integration, an adaptive grid was implemented in order 
to perform calculations with a smaller step where neces¬ 
sary. An algorithm which enabled us to make the grid 
denser in the rt-direction was implemented, as it was suf¬ 
ficient for the conducted analyses [26]. A local trunca¬ 
tion error indicator was the function where 

is a difference between the values of r for the fixed 
V at two adjacent hyper surfaces of constant u. It was 
used to mark areas where the denser mesh was required 
because it changed its value significantly in adequate re¬ 
gions [27] . 

Boundary conditions of the process were determined 
along the line u = v. The spacetime containing a spher¬ 
ically symmetric shell of matter is fiat inside the shell 
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FIG. 1. Families of the scalar field initial profiles 4>(0,i’). The amplitudes were set as equal to 0.1, 5.5, 5.0, 0.75, 147.5, 
0.115, 0.185, 700 and 0.135 for families FI to F9, respectively. 


and at large radii from it. This fact justified an assump¬ 
tion that the line u = v was not affected by the pres¬ 
ence of collapsing matter and its vicinity was nearly 
flat, which gave r = 0 along u = v. Using the equa¬ 
tion E3 the boundary condition A = 0 was obtained, 
which in combination with equations U3 and U4 led 
to / = —g = — The requirement = 0 

along u = V guaranteed the flattening of the functions 
nearby r = 0 , which made the numerical analysis pos¬ 
sible. Since r was a nonlinear function of both null co¬ 
ordinates, the three-point regressive derivative method 


with a variable step was used to implement this bound¬ 
ary condition, apart from the first point where the Euler’s 
method was used. Due to the form of the equation S 
the relation 77 = 0 held along u = v. The definition 
of 77 gave X = y along this line. Boundary conditions 
for the quantities evolving along the Tx-coordinate were 
obtained according to the equation (23) with auxiliary 
quantities (26) taken as equal to the corresponding func¬ 
tions at r = 0 . 


Initial conditions were formulated along an arbitrarily 
chosen null hypersurface, denoted as tx = 0. The assump- 
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tion of a negligible spacetime curvature in the area where 
the collapse begins justified the condition d{0,v) = 0, 
which finally fixed the remaining gauge freedom in v- 
coordinate. The flatness of spacetime in the vicinity 
of the line u = which was discussed above, gave 
a{0^v) = 1. As was mentioned in Section [lH the only 
arbitrary initial condition was the profile of the scalar 
held ^ (0, v). The quantity y (0, v) was computed analyt¬ 
ically using the equation VQ. The evolution of functions 
g, r, f and x was followed using the three-point Simp¬ 
son’s method except the first point, where the Newton’s 
method was implemented. 

The accuracy of the program was tested in two ways, 
namely the convergence of the obtained results and mass 
conservation in spacetime were checked. In order to mon¬ 
itor the convergence of the code computations for all 
the examined initial profiles (presented in Fig. were 
conducted on four grids, whose integration steps were 
multiples of 10“^. A step of a particular grid was twice 
the size of a denser one. The convergence was examined 
on rt = const, hyper surfaces chosen for each profile arbi¬ 
trarily. The values of retarded time were selected so that 
the hypersurface was situated close to the event horizon, 
but in the region where the adaptive grid was yet inactive, 
which was necessary for performing a proper comparison 
of results. 

The held functions along the selected hypersurfaces 
of constant u are shown in Fig. The areas in which 
the differences among functions obtained on various 
grids were most significant were magnified. The maxi¬ 
mum discrepancy between the finest and coarsest grids 
equal to 2.6% was observed for the profile P5. The re¬ 
maining profiles displayed even better behavior, because 
the above-mentioned differences did not exceed 1.3%. 
Figure [^presents the linear convergence of the numerical 
code. The maximal divergence between profiles obtained 
on two grids with a quotient of integration steps equal 
to 2 and their respective doubles was 3.4%. As expected, 
the errors became smaller linearly as the grid density 
increased. The convergence analysis revealed that both 
the proposed algorithm and the numerical code were ad¬ 
equate for solving the system of equations (11) - (21) 
with all investigated initial profiles PI to P9. 

The mass conservation in the obtained spacetimes 
was investigated with the use of a notion of Hawk¬ 
ing mass [28]. It is a quasi-local mass function, which 
for the Schwarzschild-type spacetimes in double null co¬ 
ordinates can be written as m 


( \ /^i , 


= %! + >) ( 27 ) 



V 



V 



V 


and it describes the mass included in a sphere of the ra¬ 
dius r (rt, v). The Hawking mass as a function of retarded 
time along the line v = 7.5, which was a maximal value 
of advanced time achieved numerically, for all the inves¬ 
tigated families of initial profiles is presented in Fig. [^ 
Since during the evolution the matter was scattered by 
the gravitational potential barrier when the collapsing 


FIG. 2. (Color online) The convergence of the scalar 
field function during the collapse initiated with the exam¬ 
ined profiles. The functions 4> (r^) were plotted for evolu¬ 
tions conducted with integration steps, which were multiples 
oi h— 10“^, along hypersurfaces of constant u equal to 0.9, 
0.1504, 0.2752, 0.3, 0.3752, 1.1, 0.1504, 1.0 and 0.7504 for pro¬ 
files PI to P9, respectively. 
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FIG. 3. (Color online) The linear convergence of the code. 
The differences between field functions calculated on grids 
with different integration steps (multiples of h = 10“^) 
and their doubles were obtained for prohles PI to P9 along 
the same hypersurfaces of constant u as in Fig. 
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FIG. 4. (Color online) Hawking masses calculated along 
V = 7.5, m/, as functions of retarded time, u, for spacetimes 
obtained with initial profiles PI to P9. 


shell approached its gravitational radius, the mass con¬ 
servation law was not satisfied in the entire domain of in¬ 
tegration. The effect of the outgoing mass flux was negli¬ 
gible except the vicinity of the event horizon. The devia¬ 
tion from mass constancy increased with advanced time, 
as the event horizon was approached. The maximal per¬ 
centage deviations from mass conservation for the initial 
profiles PI to P9 are presented in Table |lj In the majority 
of cases they were smaller than 1.0%. The most signifi¬ 
cant deviation, exceeding 16%, was observed for the pro¬ 
file P5. Table [I| also presents the values of u correspond¬ 
ing to the event horizon locations and to the percentage 
deviations not exceeding 5%. The analysis of mass con¬ 
servation in spacetime led to the conclusion that the be¬ 
havior of matter investigated numerically for all initial 
profiles was correct within the domain of integration. 
Only the results obtained for the profile P5 should be 
treated with reserve, but they were also not excluded 
from further analysis by the performed code accuracy 
checks. 


During the computations the fulfillment of equalities 
f,v = 9,u and x^y = which stem from the substi¬ 
tutions (10), was monitored for all the examined pro¬ 
files. The checks were performed in both null directions 
along V = Vf and constant whose values for the par¬ 
ticular profiles were gathered in the caption of Fig. 
The necessary derivatives were calculated with the use 
of the three-point regressive derivative method. The re¬ 
lations f^y = g^u and x^y = y^u were satisfied with 
an accuracy of 5.8% and 2.5% along the constant u 
hyper surfaces. The accuracy dX v = const, was equal 
to 2.5% and 0.1%, respectively. The discrepancies in¬ 
creased as r ^ 0 at = const, and as tx ^ ueh at con¬ 
stant V. The obtained accuracies of the constraints con¬ 
firmed the correctness of the code. 
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TABLE L The maximal percentage deviations from mass conservation, DEV, the values of retarded time corresponding 
to the event horizon location, ueh, and to the percentage deviations from mass conservation not exceeding 5%, for the in¬ 
vestigated initial profiles. 


Profile 

PI 

P2 

P3 

P4 

P5 

P6 

P7 

P8 

P9 

DEV [%] 

0.65 

0.7 

7.85 

0.62 

16.27 

0.97 

0.21 

2.53 

0.13 

UeH 

1.0609 

0.2154 

0.3421 

0.4697 

0.4954 

1.2729 

0.2858 

1.2367 

0.9331 

U5 

- 

- 

0.2631 

- 

0.1309 

- 

- 

- 

- 


IV. RESULTS 


The dynamical spacetime resulting from the scalar field 
gravitational collapse is either non-singular or it contains 
a black hole of a Schwarzschild-type, i.e., a spacelike sin¬ 
gularity situated along r = 0 surrounded by one horizon. 
The final outcome of the process depends on the value 
of the amplitude which determines the strength 
of a self-interaction of the scalar field. Since we were 
interested in the behavior of the scalar field in regions 
of high curvature, i.e., nearby the singularity, we fo¬ 
cused on amplitudes leading to the formation of a black 
hole in spacetime. They were selected for every family 
of initial profiles separately and their values are given 
in the caption of Fig. 

Penrose diagrams of the obtained dynamical space- 
times, which present r = const, lines in the (i;i4)-plane 
with 0 and (j) suppressed, are shown in Fig. The cen¬ 
tral singularity is situated along the non-linear part 
of the outermost r = const, line, which was bolded 
on the diagrams. It is spacelike, because the deriva¬ 
tive ^|^_Q is negative along the line. The singularity 
is surrounded by a single apparent horizon, marked red 
on the diagrams, which is situated along a hypersur- 
face, which becomes null in the region where v ^ oo. 
As an apparent horizon coincides with an event horizon 
in static spacetimes, it is postulated that its null charac¬ 
ter at infinity indicates the location of the event horizon 
in spacetime. Event horizons of the emerging spacetimes 
were presented as green dashed null lines on the Penrose 
diagrams. 

The location of an apparent horizon in the partic¬ 
ular spacetime allows us to divide the neighborhood 
of the central singularity into two parts. The first one cor¬ 
responds to the range of the '^’-coordinate, within which 
the apparent horizon changes its position along u. This 
region will be referred to as dynamical, because the ap¬ 
parent horizon is situated inside the event horizon, which 
is a characteristic property of dynamical spacetimes m- 
In this area the actual dynamics of the collapse takes 
place, before the spacetime settles down diS v oo. 
The remaining part of the vicinity of the singularity 
is related to the apparent horizon being situated along 
a null hypersurface. This region will be dubbed as non- 
dynamical. The null hyper surfaces of constant v which 
are borders between the dynamical and non-dynamical 
regions for spacetimes obtained with the use of profiles 


PI to P9 are pictured on the plots in Fig.j^as blue dashed 
lines. The respective values of advanced time, Vdym are 
presented on the adequate diagrams. 


The values of Vdyn were calculated on the basis of ex¬ 
amining the changes of u and '^’-coordinates along the ap¬ 
parent horizon. The changes of retarded time {5u) are 
constant due to the construction of the numerical grid 
(see Section III). The differences in advanced time {5v) 
on adjacent u = const, hyper surfaces vary from 6u 
to the values of an order of 10^6u. The values of 6v 
close to 6u are characteristic for the dynamical region, 
while the growth of Sv means that the apparent horizon 
settles along u = const, indicating the non-dynamical 
area. The selected value of Vdyn was the biggest value 
of advanced time corresponding to 6v ^ Itfidu. 


The condition which has to be fulfilled in order to treat 
the dynamical scalar field as a time variable is that 
the lines of constant values of the field are spacelike. 
The trajectories of ^ = const, lines in spacetimes emerg¬ 
ing from the collapse initiated with the use of the pro¬ 
files PI to P9 are presented in Fig. The spacetime 
regions, in which the lines are spacelike were marked 
gray on the plots and the dynamical areas nearby central 
singularities were magnified. The character of the par¬ 
ticular constant ^ line was determined on the ba¬ 
sis of the sign of the derivative calculated 

along it. The areas of spacelike and timelike hypersur- 
faces ^ = const, were characterized by positive and neg¬ 
ative values of the derivatives, respectively. 

It may be inferred from the diagram for the profile PI 
that in the vicinity of the singularity the lines of con¬ 
stant $ are spacelike, except for a limited area, in which 
the lines are timelike or null and which seems to have 
one boundary point situated at the singularity. The area 
is situated outside of the dynamical part of the vicinity 
of the singularity. Hence, nearby the central singularity 
the scalar field can be treated as a time variable. 


The interpretation of the results obtained for the re¬ 
maining initial profiles is analogous, having in mind two 
disparities. Firstly, the amount of the above-mentioned 
areas, whose boundaries have single common points with 
the singularity, depend on the family of the initial profile 
and these areas can be situated within the dynamical re¬ 
gions nearby the singularity (e.g., profiles P3, P7 and P9). 
Secondly, there can exist certain regions situated in non- 
dynamical parts of the neighborhoods of the singularities, 
in which the ^ = const, lines are non-spacelike (e.g., pro- 
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FIG. 5. (Color online) Penrose diagrams of dynamical Schwarzschild spacetimes formed during the scalar field gravitational 
collapse with initial profiles PI to P9. The borders between dynamical and non-dynamical regions are situated along null 
hypersurfaces denoted as Vdyn- 


files P4, P8 and P9). Despite the outlined differences, 
the results obtained for all the profiles enhance the com¬ 
prehensiveness of the conclusion that a scalar field can be 
treated as a time variable safely, at least in the dynami¬ 
cal part of the vicinity of the singularity of the evolving 
spacetime. 

Another requirement, which is crucial for treating 
the scalar field as a time variable is its monotonicity 
in the regions of interest. An inspection of the course 
of lines of constant 4> shown on the plots in Fig. al¬ 


lows us to draw a conclusion that in the regions, where 
the lines are spacelike, the changes are monotonic along 
constant u and v. Since the time coordinate is defined 
as t = in double null coordinates and the t-axis 
forms an angle of 45° with both u and v axes, the above 
conclusion can be extended to the coordinate t. The non¬ 
monotonic behavior of 4> = const, lines arises only when 
areas of spacelike and timelike or null characters are con¬ 
sidered, which is excluded regarding the aforementioned 
demand of the spacelike nature of constancy hypersur- 
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FIG. 6. (Color online) Lines of constant $ in spacetimes developed during the gravitational collapse of a scalar field with initial 
profiles PI to P9. Gray areas indicate spacetime regions, in which the hypersurfaces $ = const, are spacelike. The dynamical 
regions in the vicinities of central singularities were magnified. 


faces of <F. 

An important aspect of the conducted analysis, which 
stems from the fact that the numerical methods were em¬ 
ployed, is whether the singularity was approached closely 
enough in order to draw conclusions about the behavior 
of the evolving field in its close proximity. The estimation 
of the singularity closeness which was achieved during 
numerical computations is presented in Table [ll| It was 
made along the whole central singularity on the basis 
of the relation between the value of the radial function r 


which refers to the location of the event horizon in space- 
time, veh^ and its lowest value obtained in the course 
of simulations, Tmin- It turned out that the minimal val¬ 
ues of the radial function are at least six orders of mag¬ 
nitude smaller in comparison to the ones which refer 
to the event horizon location. This defines the relative 
degree of numerical closeness to the singularity obtained 
during the evolution of the scalar field. 

It should be stressed that the presented analysis is valid 
for the investigated families of profiles presented in Fig.[^ 
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TABLE IL Computational closeness to the singularity, teh denotes the value of the radial function related to the event horizon 
location in spacetime, while Tmin is the minimum value of r, which was achieved numerically in the vicinity of the central 
singularity. 


Prohle 

PI 

P2 

P3 

P4 

P5 

P6 

P7 

P8 

P9 

Teh 

0.16003 

0.23711 

0.12424 

0.51198 

0.20897 

0.36999 

0.78311 

0.41127 

0.62098 

T'min ‘ 10 

11.23 

2.60 

1.88 

0.48 

0.20 

18.72 

1.74 

0.36 

52.65 

'^min , tO'^ 

rEH 

70.17 

10.97 

15.13 

0.94 

0.93 

50.60 

2.22 

0.88 

84.79 


which means that it is independent of the selected values 
of initial amplitudes A^. Due to the widely accepted 
property of the scalar field dynamical collapse, which is 
its independence of the type of the initial profile [30] . 
the conclusions can be also extended to other profiles 
having properties described at the end of Section |n| 

V. CONCLUSIONS 

The possibility of investigating dynamics of quantum 
gravitational systems with respect to one of their matter 
components has been successfully exploited in various re¬ 
searches within quantum gravity. In our paper we made 
an attempt to justify using a scalar field as an intrin¬ 
sic ’clock’ during investigations of dynamical systems. 
For this purpose, we analyzed a dynamical evolution 
of a self-interacting scalar field within the framework 
of general relativity, for a collection of nine different ini¬ 
tial configurations of the field. 

As a result of each evolution a Schwarzschild space- 
time was obtained. It contained a spacelike singularity 
surrounded by an apparent horizon, which for all inves¬ 
tigated initial conditions was spacelike in the dynamical 
part of the spacetime. This is an interesting observation 
when associated with the notion of dynamical horizons, 
which are quasi-local spacelike horizons used to investi¬ 
gate numerous diversified aspects of physics of evolving 
black holes m- 

The interpretation of the results was concentrated 
on the issue, whether the evolving scalar field fulfills 
the necessary requirements to be regarded as a time vari¬ 
able. Our studies led to the following conclusions. 

1. The lines of constant scalar field are spacelike 
in the dynamical part of the vicinity of the cen¬ 
tral spacelike singularity formed during the dynamical 
gravitational collapse for all the investigated families 
of initial profiles. For families F4, F8 and F9 there ex¬ 


ist regions situated in non-dynamical areas, in which 
the lines are timelike or null. 

2. The scalar field values change monotonically with v 
and t coordinates in the regions, in which the hyper¬ 
surfaces of constant 4> are spacelike. 

3. The distance from the central singularity achieved nu¬ 
merically was estimated in relation to the location 
of the event horizon in spacetime. Although numeri¬ 
cal computations did not allow us to trace the scalar 
field evolution precisely up to the central singularity, 
the relative distance from it was at least six orders 
of magnitude smaller than the event horizon location, 
which means that it was highly satisfactory in the con¬ 
text of investigating the behavior of the field in its 
vicinity. 

The above findings revealed that the scalar field can 
be treated as a time variable in dynamical parts of evolv¬ 
ing spacetimes for certain, while the non-dynamical parts 
should be treated cautiously. It should be stressed 
that the presented outcomes are valid for the mass¬ 
less scalar field coupled to gravity only. In the case 
of more complicated gravitational systems, using matter 
which is coupled not only to gravity, but also to addi¬ 
tional matter fields, as a time measurer requires separate 
investigations, probably adapted to the specific prob¬ 
lem. Some indications that matter-geometry systems, 
which include non-minimally coupled Brans-Dicke field 
and a real or complex scalar field can admit intrinsic 
’clocks’ can be found in [32l[33] . 
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